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Abstract 

We examine the pion realization of the short ranged supersymmetric (SUSY) 
mechanism of neutrinoless double beta decay (O^/3/3-decay) . It originates 
from the i?-parity violating quark-lepton interactions of the SUSY extensions 
of the standard model of the electroweak interactions. We argue that pions 
are dominant SUSY mediators in Oz^/3/3-decay. The corresponding nuclear 
matrix elements for potentially O^/3/3-decaying isotopes are calculated within 
the proton-neutron renormalized quasiparticle random phase approximation 
(pn-RQRPA). We define those isotopes which are most sensitive to the SUSY 
signal and outlook the present experimental situation with the O^/3/3-decay 
searches for the SUSY. Upper limits on the i?-parity violating 1st generation 
Yukawa coupling A' m are derived from various Oz^/3/3-experiments. 

12.60.Jv, 11.30.Er, 11.30.Fs, 23.40.Bw 
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I. INTRODUCTION 



The observation of neutrinoless nuclear double beta decay 

(A, Z) — >(A,Z + 2) + 2e~ (1) 

(Ovflfi) would undoubtedly indicate the presence of the new physics beyond the standard 
model (SM) of electroweak interactions. However, as yet there is no experimental evidence 
for this lepton-number violating (AL = 2)exotic process. On the other hand non-observation 
of Oz//3/3-decay at certain experimental sensitivity allows one to set limits on some parameters 
of the new physics. An unprecedented accuracy and precision of the modern 0z//3/3-decay 
experiments allows one in certain cases to push these limits far out of reach of the accelerator 
and the other non-accelerator experiments. 

A well known example is given by the upper limit on the light effective Majorana neutrino 
mass (m v M ). From the 0i//9/3-decay experiments [0] it was found (m u M ) < 0(1.1 eV) 0. 
Recall that the Majorana neutrino mass term violates the lepton number AL = 2. This 
is exactly what is necessary for (V/?/3-decay to proceed via the virtual neutrino exchange 
between the two neutrons. In this case the 0z//?/3-decay amplitude is proportional to (m^). 

The Majorana neutrino exchange is not the only possible mechanism of Oz//5/5-decay. 
The lepton-number violating quark-lepton interactions of the R-parity non-conserving su- 
persymmetric extensions of the SM SUSY) can also induce this process ||- |J. .R-parity 
is a discrete multiplicative Z 2 symmetry defined as R p = (_i) 3B + i + 2S ' ) where S, B and L 
are the spin, the baryon and the lepton quantum number. At the level of renormalizable 
operators R-parity can be explicitly violated by the trilinear and the bilinear terms in the 
superpotential and in the soft SUSY breaking sector. 

The impact of the R-parity violation on the low energy phenomenology is twofold. First, 
it leads the lepton number and lepton flavor violating interactions directly from the trilin- 
ear terms. Second, bilinear terms generate the non-zero vacuum expectation value for the 
sneutrino fields (z/j) ^ and cause neutrino-neutralino as well as electron-chargino mixing. 
The mixing brings in the new lepton number and lepton flavor violating interactions in the 
physical mass eigenstate basis. 

The implications of the trilinear and the bilinear terms for the Oi//3/3-decay were previ- 
ously considered in Refs. ||- and in Ref. respectively. The Oz//3/3-decay proved to be 
very sensitive probe of the new interactions predicted in the Jfi p SUSY ||- ||. 

In this paper we return to the phenomenology of the trilinear terms and perform a 
comprehensive analysis of their contribution to the 0z//9/9-decay, paying special attention 
to the hadronization of the corresponding quark interactions and to the nuclear structure 
calculations. We will show that for the case of the trilinear terms the stage of hadronization 
plays especially important role in derivation of the short-ranged Tjl v SUSY mechanism of 
0^/?/9-decay. 

In our paper Ref. |]] we had considered the two-pion realization of the underlying AL = 2 
quark- level Oi//3/3-transition dd — > uu+2e~. It was found that the corresponding contribution 
to O^/3/3-decay absolutely dominates over the conventional two nucleon mode realization. In 
this paper we generalize the previous treatment of the hadronization of the AL = 2 quark 
operators. 
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Searching for tiny effects of the physics beyond the SM in 0uf3f3-decay requires a reli- 
able treatment of the nuclear structure as well. In this paper we present the results of our 
calculations within the proton-neutron renormalized Quasiparticle Random Phase Approx- 



imation (pn-RQRPA) H, [ lQ[ . We are listing the nuclear matrix elements for the lfl p SUSY 
mechanism of the 0i//?/9-decay and their specific values for the experimentally interesting 
isotopes. 

We introduce new characteristic of 0z//3/3-decaying isotope: its sensitivity to the lflp 
SUSY signal. This characteristic depends only on the corresponding nuclear matrix ele- 
ment and the kinematical phase-space factor. We calculate the SUSY sensitivities of all 
experimentally interesting isotopes and on this basis estimate prospects for SUSY searches 
in 0z//?/5-experiments. Applying our approach, we determine presently most successful 
0i//3/3-experiment which establishes the most stringent constraints on the R-parity violating 
Yukawa couplings. 

The paper is organized as follows. In the next section we shortly outline the minimal su- 
persymmetric standard model with the explicit R-parity violation (]ji p MSSM) and show the 
effective AL = 2 Lagrangian which describes in this model the quark-level Ozz/^-transition. 
In section 3 we derive the corresponding effective Lagrangian at the hadronic level in terms 
of the meson and the nucleon fields. Section 4 is devoted to derivation of the nuclear 0vf3f3- 
transition operators in non-relativistic impulse approximation and to calculation of their 
matrix elements in the pn-RQRPA approach. Section 5 deals with the constraints on ]ji p 
Yukawa couplings from various Oz//?/?-experiments. 



II. $ P SUSY INDUCED AL = 2 QUARK-LEPTON INTERACTIONS. 



Let us shortly outline the minimal supersymmetric standard model with the explicit 
R-parity violation (R p MSSM). 

For the minimal MSSM field contents the most general gauge invariant form of the 
renormalizable superpotential is 



w = w R + W J 



v 



where the R p conserving part has the standard MSSM form [[11] 

W Rv = h L H x LE c + h D HiQD c + h v H 2 QU c + ^H X H 2 . 



(2) 



(3) 



Here L, Q stand for lepton and quark doublet left-handed superfields while E c , U c , D c for 
lepton and up, down quark singlet superfields; H 1 and H 2 are the Higgs doublet superfields 
with weak hypercharges Y = —1, +1, respectively. Summation over the generations is 
implied. 

The R p violating part of the superpotential (0) can be written as [12], [|13| 



= XijkULjEl + X' ijk LiQjDl + fijLjH 2 + \" jk UiDjDl, 



(4) 



The coupling constants A (A") are antisymmetric in the first (last) two indices. 

The soft supersymmetry breaking part of the scalar potential also contains the ^? p -terms 
of the form: 
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V;° p ft = KikULiEi + KikLiQA + KiiPtDpt + (5) 

+fl 2 2j L j H 2 + fiijLjHi + H.c. 

In Eqs. (fD, (H) the trilinear terms proportional to A, A', A, A' and the bilinear terms violate 
the lepton number while the trilinear terms proportional to A", A" violate baryon number 
conservation. 

It is well known that the simultaneous presence of lepton and baryon number violating 
terms in Eqs. (|), (|) (unless the couplings are very small) leads to unsuppressed proton 
decay. Therefore, either only the lepton or the baryon number violating couplings can be 
present. Certain discrete symmetries such as the B-parity [OJ, |14| may originate from the 



underlying high energy scale theory and forbid dangerous combinations of these couplings. 
Henceforth we simply set A" = A" = 0. 

The remaining R-parity conserving part of the soft SUSY breaking sector includes the 
scalar field interactions 



V R° - E m 2 Mi\ 2 + h L A L H l LE c + h D A D H l QD c + (6) 

i=scalars 

+h u A u H 2 QU c + nBH x H 2 + H.c. 

and the "soft" gaugino mass terms 



1 r 



C-GM — — - 



M X BB + M 2 W k W k + M 3 g a g a - H.c. (7) 



As usual, M 3j2 ,i denote the "soft" masses of the SU(3) x SU(2) x U(l) gauginos g,W,B 
while rrti stand for the masses of the scalar fields. The gluino g soft mass M 3 coincides in 
this framework with its physical mass denoted hereafter as m~ g = M 3 . 

As mentioned in the introduction, we concentrate on the phenomenology of the trilin- 
ear R-parity violating terms LQD C in the superpotential (^) and perform a comprehensive 
analysis of their contribution to the Oz//5/3-decay. 

These terms lead to the following AL = 1 lepton-quark operators 

£a = Kjk [f> iL d h Pl^ + d. L d k P L v. + d kR d j P R v c . - e iL d k P L u 3 (8) 
- u 3L d k P Lei - d kR u ] P R e':] + H.c. 

Here, as usual Pl,r = (1 =F 7s)/2. 

Starting from this fundamental Lagrangian, one can derive the low-energy effective La- 
grangian 0, which describes the quark-level 0^/9/9-transition dd —>■ uu + 2e~. Integrating 
out the heavy degrees of freedom, we come up with the formula: 



r 2 r 1 

V JPS J PS — J T J Tfiu 



C qe =^e(l + 75 )e' 



(9) 



These AL e = 2 lepton-number violating effective interactions are induced by heavy SUSY 
particles exchange. An example of the Feynman diagram contributing to C qe is given in Fig. 
1. A complete list of the relevant diagrams can be found in M. Compared to Ref. we 
have properly taken into account in the Lagrangian (|9]) the contribution of the color octet 
currents. 
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The color-singlet hadronic currents in Eq. @ are 



Jps = Jp + Js, Jp = u a ^d ai J s = u a d a , = u a a^{l + j 5 )d 



(10) 



where a is the color index and a^ u = (i/2)[-f tJ- , j"]. 

The effective lepton-number violating parameters rj PS , r] T depend on the fundamental 
parameters of the MSSM and can be written in the form 



PS 



V X e +V x f + Vx + Vs + H 



where we denoted 

Vg 
Vx 
Vxe 
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In Eqs. (p~3|)- (P~7|) we used the standard notations a 2 = g 2 /(4:ii) and a s = gf/(47r) for 
the SU{2)l and SU(3) C gauge coupling constants. We also denoted the gluino g and the 
neutralinos Xi masses as rrig and m Xi respectively. The Majorana neutralinos Xi are linear 
combinations of the gaugino and higgsino fields 



Xi 



M iX B + M l2 W 3 + NaHl + NaHl 



(18) 



Here W 3 and B are the neutral SU(2)l and U(l) gauginos while H®, are the higgsi- 
nos which are superpartners of the two neutral Higgs boson fields H® and H® with weak 
hypercharges Y — — 1, +1, respectively. 

The matrix Mij, introduced in Eq. (|T8"D, rotates the 4x4 neutralino mass matrix M x to 



the diagonal form Diag[m Xi }. We define these matrix as in Ref. []lT 
Neutralino couplings are defined as [|TT|] 

e w (V0 = -T 3 (il>)Af i2 + tan0 w (T 3 (^) - Q(i/>))N a , 



(19) 
(20) 



Here Q and T 3 are the electric charge and the weak isospin of the fields ip = u,d, e. 
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III. HADRONIZATION OF $ P SUSY QUARK-LEPTON INTERACTIONS 



The next step deals with reformulation of the quark-lepton interactions in Eq. (^) in 
terms of the effective hadron-lepton interactions. This is necessary for the subsequent nuclear 
structure calculations. 

There are the two possibilities of hadronization of the effective Lagrangian C qe in Eq. 
One can place the four quark fields in the two initial neutrons and two final protons 
separately. This is the conventional 2N-mode of 0z//?/3-decay shown in Fig. 2(a). Then 
nn — > pp + 2e~ transition is directly induced by the underlying quark subprocess 

dd->im + 2e~. (21) 

In this case the nucleon transition is mediated by the exchange of a heavy supersymmetric 
particle like the gluino g with the mass rrig > lOOGeV. Therefore, the two decaying neutrons 
are required to come up very closely to each other what is suppressed by the nucleon-nucleon 
short range repulsion. 

Another possibility is to incorporate quarks involved in the underlying Jji v SUSY transi- 
tion in Eq. ([H|) not into nucleons but into two virtual pions || or into one pion as well as 
into one initial neutron and one final proton. Now the nn —> pp + 2e~ transition is mediated 
by the charged pion-exchange between the decaying neutrons, as shown in Fig. 2(b,c). This 
is what we call the one- and two-pion modes of Ov/3/3-decay. Since the interaction region 
extends to the distances ~ l/m n this mode is not suppressed by the short range nucleon- 
nucleon repulsion. An additional enhancement of the 7r-modes comes from the hadronization 
of the Ift p SUSY quark-lepton vertex operator in Eq. ([]) as discussed below. In Ref. it was 
shown that the two-pion mode absolutely dominates over the 2N-mode. In what follows, we 
are arguing that it dominates over the one-pion mode as well. 

The effective hadronic Lagrangian taking into account both the nucleon (p, n) and tt- 
meson degrees of freedoms in a nucleus can be written as follows: 

G 2 

Che = C 2N + C 2n + £ l7r + C s = — F - pT®n ■ pT®n • e(l + 75 )e c - (22) 

2m„ 



r 2 

mt 



2m 

p 



mla 27T (-7T ) - a l7T p ij 5 n ■ n~ 



e(l + 7 5 )e c + g s p ij 5 n vr" 



Here C s stays for the standard pion-nucleon interaction with the coupling g s = 13.4 ± 1 
known from experiment. The lepton-number violating terms £ 2 n, £in, £-2n generate the 
conventional two-nucleon mode, the one and two pion-exchange modes of the 0z//?/3-decay 
respectively. The corresponding diagrams are presented in Fig. 2. 

The two-nucleon mode term L 2 m with different operator structures had been consid- 
ered in HH within the _^ p MSSM. As was already mentioned this term gives the sub-dominant 
contribution to O^/3/3-decay in comparison with the pion terms @. Therefore, in this paper 
we concentrate on the effect of the pion-exchange contribution generated by the terms C\ n 



and C 2w in Eq. (H). 

The basic parameters a 2n and a\ n of the Lagrangian Che m Eq- ([22] ) can be approximately 
related to the parameters of the quark-lepton Lagrangian C qe in Eq. (|9|), using the on- mass- 
shell "matching conditions" [El 
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< 7r + ,2e-|£g e |7r~ > = < -K + ,2e'\C 27T \^ >, (23) 
< ti + ,p,2e~\C qe \n > = < ii + ,p,2e~\Ci n \n > . (24) 

In order to solve these equations we apply the widely used factorization and vacuum domi- 
nance approximations [15| for the matrix elements of the products of the two quark currents. 
Then we obtain, taking properly into account the combinatorial and color factors: 

(n + \J PS J PS \n-) « ^(n + \J P \0)(0\J P \n-), 

(tt+IJ^J^Itt-) « -4(7r+|Jp|0)(0|Jp|7r->, (25) 

(p\JpsJps\nir-) « -(p|J F |n)(0|J P |7T-), (26) 

(plJ^Jr^InO « -4{p|J P |7i)(0|Jp|7r-). (27) 

Here we applied the equalities 

< o| J s \n >=< 0|J| iy |7r(p 7r ) >= 0. (28) 

The scalar matrix element vanishes due to the parity arguments, the tensor one vanishes 
due to J^ v = —J^ 1 and impossibility of constructing an antisymmetric object having only 
one external 4- vector p n . 

We also use the following relationships for the hadronic matrix elements: 

2 

(0|« T5 d|7r-) = iV2~U m * = imlK, (29) 
m u + m d 

(p\u^ 5 d\n) = F P {p\p-f b n\n) . (30) 
where f n = 0.668 m n . For the nucleon pseudoscalar constant Fp we take its bag model value 



4.41 from Ref. 16 



In this approximation we solve the "matching conditions" in Eqs. (|23f)-(|24|) and determine 
the coefficients in Eq. ( f22"|) 

OfcTT = C kn ^[l] T + ^7] PS ], (31) 

o o 

with 

cm = h w F P « 132.4, c 2n = ~hl « 170.3. (32) 

Here we accepted the conventional values of the current quark masses m u = 4.2 MeV, 
= 7.4 MeV. The large ratio of the pion mass to the small current quark masses provide 
an additional enhancement factor of the pion mechanism as mentioned before the Eq. (p2| ) . 
Thus, we have obtained the approximate hadronic "image" £h e of the fundamental quark- 
lepton Lagrangian C qe given in Eq. @. 
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IV. NUCLEAR MATRIX ELEMENTS OF $ P SUSY INDUCED 

Oi//3/3-TRANSITION 



Starting from the Lagrangian Che m Eq. (p2j) it is straightforward to calculate the Qv(3(3- 
nuclear matrix element 

< (A,Z + 2),2e~\S - 1\{A,Z) > =< (A, Z + 2), 2e~\Texp[i J d A xC he (x)]\(A, Z) > (33) 

The nuclear structure is involved via the initial (A,Z) and the final (A, Z+2) nuclear states 
having the same atomic weight A, but different electric charges Z and Z+2. The nucleon- 
level diagrams, which correspond to the leading order contributions to the amplitude in Eq. 
([33]), are given in Fig. 2. 

The standard framework for the calculation of this nuclear matrix element is the non- 
relativistic impulse approximation (NRIA) |F7| . 

The final result for the half-life of Oz//3/3-decay in + — > + channel with the two outgoing 
electrons in the S-state, regarding all the three above-described possibilities of hadronization, 
reads 

[T 1/2 (0 + ^0 + )]- 1 = (34) 



01 



T 

V 



■ Mf + ( V PS - V T ) ■ Mf + 1(V T + lv PS ) (^M 1 - + M 2 -) 



Here Goi is the standard phase space factor tabulated for various nuclei in Ref. |18| . The 
nuclear matrix elements governing the sub-dominant two-nucleon mode were presented 
in Ref. 0. As was already mentioned its contribution can be safely neglected. The one- 
and the two-pion modes nuclear matrix elements M ln and M 2lT we write down in the form 



'^\ 2 ^ a ^( M k J T + M^) (35) 



Here, = 850MeV is the mass scale of the nucleon form factor. 

The structure coefficients a kn in Eq. ( |35l) are related to the coefficients introduced in 
Eq. ( |32"D and have the following explicit form 

a l7T = -6g 3 KpF P w -4.4 • 10~ 2 , (36) 
«^ = ^> 2 p^2.0-10- 1 , (37) 



where 



with Ja ~ 1.261 being the axial nucleon coupling. 

The two types of the Gamow- Teller and tensor nuclear matrix elements are given by the 
expressions 



S 



M k J T = 


(°/ 


E 


= 


(0/ 


E 



,(k), 



i? 



(39) 
(40) 



where 



rj is a coordinate of the ith nucleoli and R = tqA 1 ^ 3 is the mean nuclear radius, with 
r = 1.1 fm. 

The structure functions Fi(x n ) and F2 (x n ) have their origin in the integration of the 
pi-meson propagators and take the following form: 



F x (1) (x) = e~ x , F 2 (1) (x) = (3 + 3x + x 2 )^, (42) 



e 



x 



F{ 2 \x) = {x- 2)e~*, F 2 (2) (x) = (x + l)e~ x . (43) 

We calculate the nuclear matrix elements within the proton-neutron renormalized Quasi- 
particle Random Phase Approximation (pn-RQRPA) PJIOfl. This nuclear structure method 
has been developed from the proton-neutron QRPA (pn-QRPA) approach, which has been 
frequently used in the Ou/3/3-decaj calculations. The pn-RQRPA is an extension of the pn- 
QRPA by incorporating the Pauli exclusion principle for the fermion pairs. The limitation of 
the conventional pn-QRPA is traced to the quasiboson approximation (QBA), which violates 
the Pauli exclusion principle. In the QBA one neglects the terms coming from the commu- 
tator of the two bifermion operators by replacing the exact expression for this commutator 
with its expectation value in the uncorrelated BCS ground state. In this way the QBA 
implies the two-quasiparticle operator to be a boson operator. The QBA leads to too strong 
ground state correlations with increasing strength of the residual interaction in the particle- 
particle channel what affects the calculated nuclear matrix elements severely. To overcome 
this problem the Pauli exclusion principle has to be incorporated into the formalism ||, 
TOfl in order to limit the number of quasiparticle pairs in the correlated ground state. The 



commutator is not anymore boson like, but obtains corrections to its bosonic behavior due 
to the fermionic constituents. The pn-RQRPA goes beyond the QBA. The Pauli effect of 
fermion pairs is included in the pn-RQRPA via the renormalized QBA (RQBA) [TO, i.e. 
by calculating the commutator of two bifermion operators in the correlated QRPA ground 
state. The RQBA was applied to the 2z//?/3-decay in Ref. || and to the Oz//3/3-decay for the 
first time in Ref. f]T0|| . Now it is widely recognized that the QBA is a poor approximation 
and that the pn-RQRPA offers the advantages over pn-QRPA. Let us stress that there is 
no collapse of the pn-RQRPA solution for a physical value of the nuclear force and that 
the nuclear matrix elements have been found significantly less sensitive to the increasing 
strength of the particle-particle interaction in comparison with QRPA results [f2] . Thus, the 
pn-RQRPA provides significantly more reliable treatment of the nuclear many-body problem 
for the description of the Qv(3(3 decay. 

For numerical treatment of the 0z//9/9-decay matrix elements given in Eqs. (|3^) and ([40]) 
within the pn-RQRPA we transform them by using the second quantization formalism to 
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the form containing the two-body matrix elements in the relative coordinate. One obtains 



pnp' n' 



< o; || [cp n ,]j || J^m/ >< rm f \ r mi >< r mi \\ [ c +g n ]j || of > . (44) 

Oij represents the coordinate and spin dependent part of the two body transition operators of 
the 0i//9/9-decay nuclear matrix elements in Eqs. ( |3"§D and fl4"0j). The short-range correlations 
between the two interacting nucleons are taken into account by the correlation function 

f{r) = l-e~ ar \l-br 2 ) with a = 1.1 fm 2 and 6 = 0.68 fm 2 . (45) 

The one-body transition densities and the other details of the nuclear structure model are 
given in 

The calculated nuclear matrix elements for the Oz//3/3-decay of various isotopes within 
the pn-RQRPA are presented in Table 1. The considered single-particle model spaces both 
for protons and neutrons have been as follows: i) For A=76, 82 the model space consists of 
the full 2 — 4hu major oscillator shells, ii) For A=96, 100, 116 we added to the previous 
model space I/5/2, I/7/2, 0/t 9 / 2 and 0hu/ 2 levels, iii) For A=128, 130, 136 the model space 
comprises the full 2 — 5hu major shells, iv) For A=150 the model space extends over the 
full 2 — 5hu shells plus Oin/2 an d 0zi 3 / 2 levels. 

The single particle energies were obtained by using Coulomb corrected Woods Saxon 
potential. The interaction employed was the Brueckner G-matrix which is a solution of the 
Bethe-Goldstone equation with the Bonn one-boson exchange potential. Since the considered 
model space is finite, we have renormalized pairing interactions by the strength parameters 
dpp and d nn [[0| to the empirical gaps defined by Moeller and Nix [p0|1 . The particle-particle 



and particle-hole channels of the G-matrix interaction of nuclear Hamiltonian H have been 
renormalized with parameters g pp and g p h, respectively. The nuclear matrix elements listed 
in the Table 1 have been obtained for the g p h = 0.80 and g pp = 1.0 

The following note is in order. According to our numerical analysis, variations of the 
nuclear matrix elements presented in Table 1 do not exceed 20% within the physical region 
of the nuclear structure parameter g pp (0.8 < g pp < 1.2). 

As seen from the Table 1 M 2n is significantly larger than M l7T . It is partially due to the 
mutual cancellation of the M}? T and Mf 71 " in the construction of M ln in Eq. (^) (see Table 
1) and due to the suppression of the structure coefficient a ln in comparison with a 2n . Thus, 
we conclude that the two-pion mode contribution to 0^/?/?-decay (Fig. 2(c)) dominates both 
over the one-pion (Fig. 2(b)) and the two-nucleon contributions (Fig. 2(a)). The dominance 
of the two-pion mode over the two-nucleon one was previously proven in Ref. 0. 



V. CONSTRAINTS ON $ P SUSY FROM Oi//3/3-EXPERIMENTS. COMPARATIVE 

ANALYSIS. 

Having all the quantities in the 0i//3/9-decay half-life formula Q3~4]) specified we are ready 
to extract the limits on the Tjl p parameters from non-observation of the Oz//5/5-decay. 
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The experimental lower bound T x m for the half-life of a certain isotope Y provides the 
following constraint on the effective $ p SUSY parameters 



+ 3 V 



PS\ 



10" 



10 24 years 

rpexp > 

J l/2 



(46) 



Here we introduced the SUSY sensitivity COO °f a 0z//9/?-decaying isotope Y 



C(Y) = 10 5 \-M ln + M 



r2n\ 



a 



01- 



(47) 



The quantity COO is an intrinsic characteristic of an isotope Y depending only on the nuclear 
matrix elements M ln , M 2n and on the phase space factor Gq\. The large numerical values 
of the SUSY sensitivity C defined in ([47]) correspond to those isotopes within the group 
of /?/9-decaying nuclei which are the most promising candidates for searching SUSY in the 
O^/3/3-decay. 

The numerical values of C(^0 calculated in the pn-RQRPA are presented in the Table 1 
and displayed in Fig. 3 in the form of a histogram. It is seen that the most sensitive isotope 
is 150 Nd, then follows 100 Mo. 

It is understood that the SUSY sensitivity C cannot be the only criterion for selecting 
an isotope for the 0^/3/3-experiment. Other microscopic and macroscopic properties of the 
isotope are also important for building a Oz//3/5-detector. 

The current experimental situation in terms of the accessible half-life and the correspond- 
ing upper limit on the effective SUSY parameter r] SUSY is presented in Table 2. We conclude 
that the best upper limit on the Ifip SUSY parameter i] SUSY has been established by the 
Heidelberg- Moscow experiment |]J . We denote this limit as Tj^f (H — M) . For comparison 
in the bottom row of the Table 2 we show the lower half-life limits T^(rjg~J^), which must 
be reached by 0i//3/3-experiments with the other nuclei to reach this presently best constraint 
Vsusy — ^Tusy ~ M) on the lfi v SUSY. The result of this comparison is illustrated in Fig. 
3. 

Using the values of ff^ SY from the Table 2 one can easily calculate the corresponding 
constraints on the A' m parameter. There are two types of the constraints for each value of 

Vg^SY P arame t er 



(- 



m?. 



/I 



1/2 



A'm < 12.5 



exp 

SUSY VlOOGeW VlOOGeW ' 

2 / /m \ 1/2 



^exp 



nir. 



SUSY 



X(- 



VIOOGeV/ VlOOGeV; 



(48) 
(49) 



These formulas are derived from Eqs. QllD- (|17D ) applying a widely used ansatz of the uni- 



versal squark q mass rrig at the weak scale m s 



m d 



ni, 



This approximation is well 



motivated by the constraints from the flavor changing neutral currents. Formula ([49]) takes 
into account only the lightest neutralino contribution with the mass m x . We also assume 
absence of spurious compensations between terms of different nature such as the g — q — q and 
X — e — e. The running QCD coupling constant a s (Q) has been taken at the scale Q = lGeV 
with the normalization on the world average value a s (Mz) = 0.120 | 29| . The second limit 
in Eq. (|49|) is derived with the additional assumptions that the lightest neutralino is B-ino 
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dominant and that rriq > mg/2. Both these assumptions are phenomenologically reasonable, 
although they must not be always correct. 

The best constraint from the Heidelberg-Moscow experiment ill] is 

These limits are very strong and, as it was already pointed out in Ref. ||- |7|], lie beyond 
the reach of the near future accelerator experiments (though, accelerator experiments are 
potentially sensitive to the other couplings than A' m ). 

To constrain the size of A' m itself one needs additional assumptions on the masses of the 
SUSY-partners. If the values of these masses would be around their present experimental 
lower limits ~ lOOGeV [29(1, one could constrain the coupling to 

A' m < 1.3 • 10" 4 . (52) 

A conservative bound can be set by assuming all the SUSY-masses being at the "SUSY- 
naturalness" bound of ~ ITeV, leading to 

A' m < 0.04. (53) 

This completes our analysis. The other details concerning the experimental prospects for 
searching for lfi p SUSY in 0z//9/?-experiments can be inferred directly from Table 2 and Fig. 
3. 



VI. CONCLUSION 

In summary, we have analyzed the general case of the pion realization for the short- 
ranged lfip SUSY mechanism of O^/5/3-decay taking into account both the one-pion and 
two-pion modes. We have shown that the two-pion mode Ifi p SUSY contribution to Qv(3(3- 
decay dominates over the one-pion mode contribution. Previously || we had proven that 
the two-pion mode dominates over the conventional two-nucleon one. We also pointed out 
that non-observation of 0z//?/9-decay casts severe limitations on the lji p SUSY extensions of 
the standard model of electroweak interaction. Although a complicated nuclear many-body 
problem needs to be solved the limits are so stringent that they overcome the uncertainties 
in the nuclear and hadronic matrix elements, leading to limits that are much stronger than 
those from accelerator and the other non-accelerator experiments. We gave the list of nuclear 
matrix elements for all the experimentally interesting isotopes and presented the so called 
SUSY sensitivities of these isotopes. These characteristic might be helpful for planning 
future searches for SUSY in 0^/3/3-decay. On this basis we compared the present status of 
the various 0^/?/9-experiments and their abilities to detect the SUSY signal. 
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TABLES 



TABLE I. Nuclear matrix elements for the pion-exchange R-parity violating SUSY mode of 
Oi^/5/5-decay for the experimentally most interesting isotopes calculated within the renormalized 
pn-QRPA. Goi is the integrated kinematical factors for + — ► + transition [18]. C(X) denotes 
according to Eq. ( [47| ) the sensitivity of a given nucleus Y to the SUSY signal. 



Nucleus 


M% T 


M% 


M 2 J T 


M$ 


M 1 * 


M 27T 


xl0 15 y 


C(T) 


76 Ge 


1.30 


-1.02 


-1.34 


-0.65 


-18.2 


-601 


7.93 


5.5 


82 Se 


1.23 


-0.87 


-1.26 


-0.57 


-23.9 


-551 


35.2 


10.8 


96 Zr 


0.77 


-1.11 


-0.85 


-0.67 


22.1 


-458 


73.6 


11.8 


W0 Mo 


1.43 


-1.73 


-1.52 


-1.05 


19.4 


-776 


57.3 


18.1 


ne Cd 


0.92 


-0.78 


-0.94 


-0.47 


-9.3 


-423 


62.3 


10.8 


128 Te 


1.25 


-1.57 


-1.40 


-0.99 


21.4 


-720 


2.21 


3.3 


i30 Te 


1.10 


-1.48 


-1.26 


-0.93 


25.1 


-660 


55.4 


14.9 


136 Xe 


0.61 


-0.84 


-0.74 


-0.54 


15.5 


-387 


59.1 


9.0 


150 Nd 


1.85 


-2.70 


-2.07 


-1.68 


56.4 


-1129 


269. 


55.6 



TABLE II. The present state of the $ p SUSY searches in /3/3-decay experiments. Ty% (present) 
is the best presently available lower limit on the half-life of the 0^/3/?-decay for a given isotope. 
rfsusY * s ^ ne corresponding upper limit on the $ p SUSY parameter. T^(wfj~Jy) is the calculated 



half-life of Oi^/5/5-decay assuming 


n — jjH—M 

'ISUSY ''SUSY 


with rjf US Y being 


• the best limit deduced from the 


Heidelberg-Moscow 


76 Ge experiment [1]. 








Nucleus 


76 Ge 


82 Se 


96 Zr 


100 Mo 


116 Cd 


Ref. 


ng 


m 


© 




m 


Ty£ (present) 


1.1 x 10 25 


2.7 x 10 22 


3.9 x 10 19 


5.2 x 10 22 


2.9 x 10 22 


> / •' / ' 

•SUSY 


5.5 x 10~ 9 


5.6 x 10~ 8 


1.4 x 10~ 6 


2.4 x 10~ 8 


5.4 x 10~ 8 


rpexp, H-M\ 
± l/2 V'SUSY > 


1.1 x 10 25 


2.9 x 10 24 


2.4 x 10 24 


1.0 x 10 24 


2.9 x 10 24 



Nucleus 



128 



Te 



130 



Te 



136 



Xe 



150 



Nd 



m 



m 



Ty£ (present) 



n exp 



rpexp ( H-M\ 
1/2 V'SUSY > 



7.7 x 10 24 
1.1 x 10~ 8 
3.1 x 10 25 



8.2 x 10 21 

7.4 x 10~ 8 

1.5 x 10 24 



4.2 x 10 23 
1.7 x 10~ 8 
4.1 x 10 24 



1.2 x 10 21 
5.2 x 10~ 8 
1.1 x 10 23 
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Figure Captions 



Fig. 1: An example of the supersymmetric contribution to Oz//3/3-decay with the gluino g 
and two squarks u in the intermediate state. 

Fig. 2: The hadronic-level diagrams for the short-ranged SUSY mechanism of (V/3/3-decay. 
(a) the conventional two-nucleon mode, (b) the one-pion exchange mode, (c) the two-pion 
exchange mode. 

Fig. 3: The SUSY sensitivity COO f° r the experimentally interesting nuclei (on the left). 
This histogram displays the corresponding numerical values in the Table 1. The histogram 
on the right illustrates the Table 2. The best presently available lower limits on the Ovftfl- 
decay half-life are denoted by the black bars. The open bars indicate the half-life limits 
TyliVsusv) t° be reached by a given experiment to reach the presently best limit on the $ p 
SUSY parameter V^Jsy established by the 76 Ge experiment [0. 
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